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Several methods for analyzing the particle spectra from extremely large solar proton 
events, called Ground-Level Enhancements (GLEs), have been developed and utilized by the 
scientific community to describe the solar proton energy spectra and have been further 
applied to ascertain the radiation exposures to humans and radio-sensitive systems, namely 
electronics. In this paper 12 GLEs dating back to 1956 are discussed, and the three methods 
for describing the solar proton energy spectra are reviewed. The three spectral fitting 
methodologies are EXP [an exponential in proton rigidity (R)[, WEIB |Weibull fit: an 
exponential in proton energy], and the Band function (BAND) [a double power law in proton 
rigidity]. The EXP and WEIB methods use low energy (MeV) GLE solar proton data and 
make extrapolations out to ~1 GeV. On the other hand, the BAND method utilizes low- and 
medium-energy satellite solar proton data combined with high-energy solar proton data 
deduced from high-latitude neutron monitoring stations. Thus, the BAND method 
completely describes the entire proton energy spectrum based on actual solar proton 
observations out to ~10 GeV. Using the differential spectra produced from each of the 12 
selected GLEs for each of the three methods, radiation exposures are presented and 
discussed in detail. These radiation exposures are then compared with the current 30-day 
and annual crew exposure limits and the radiation effects to electronics. 
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I. Introduction 

S olar protons events (SPEs) represent the single-most significant source of acute radiation exposure to humans 
and space systems during lunar and deep-space missions. Occasionally, extremely large SPEs occur, called 
Ground Level Enhancements (GLEs), that contain protons having energies extending into the GeV range. Early in 
the “Space Age era,” particle detectors on satellites had limited energy resolution capability. Malitson and Webber 1 
suggested an exponential in proton rigidity fit as a means of describing the energy spectrum of an SPE based on two 
data points: integral fluence <I> > 30 MeV and > 100 MeV (protons/cm 2 ). Beyond 100 MeV the particle energy 
spectrum was extrapolated out to several GeV. This methodology was used by the scientific community for several 
decades. Later, an energy spectrum fitting method, called a Weibull fit, was introduced by Xapsos, et al. 2 . Most 
recently, Tylka and Dietrich 3 and Atwell, et al. 4 applied a Band function 5 fit, which is a double power law in proton 
rigidity. Their method is based on actual proton data observed from medium (10-100 MeV) and high (several 
hundred MeV) proton energies (GOES MEPAD and HEP AD and other satellites) and from cascade secondary 
neutrons produced during nuclear collisions of the incoming solar protons with the Earth’s atmosphere (> 1 GV). 
These neutrons are measured at the high latitudes by a network of neutron monitoring stations. These authors have 
successfully applied the Band function fitting method to ~70 GLEs dating back to 1956. Since this method 
completely describes the entire proton energy spectrum for an event, this is the method that is strongly suggested 
that the scientific community use for future GLE SPE assessment and analysis. 

In this paper we have selected 12 GLEs from the original ~70 GLEs and present the integral and differential 
proton energy spectra ranging from 10 MeV to 5 GeV. In addition, we have taken the differential spectra and 
computed “skin” and “Blood-Forming Organ (BFO)” radiation exposures using the NASA Langley Research Center 
HZETRN 2010 6 high energy particle transport/dose code. These skin and BFO exposures for the 12 GLEs are 
compared with the current NASA 30-day and annual crew exposure limits. 

II. The Three Fitting Methods 

A. The Exponential in Proton Rigidity Method (EXP) 

The exponential in particle rigidity fitting method of Malitson and Webber 1 is based on two proton integral data 
points: >30 MeV and >100 MeV. The mathematical expression for the integral energy spectrum is: 

d> (>E) = N 0 * exp (-R/R 0 ) 

where <!> (>E) is the integral energy fluence in protons/cm 2 , N 0 is a normalization constant, R is the proton rigidity 
(proton momentum) in MV (million volts), and R 0 is the characteristic rigidity in MV. The proton rigidity R is 
related to the proton energy (MeV) by: 


R (MV) = V (E 2 - 2*m 0 *E) 

where E is the proton energy (MeV) and m Q is the rest mass of the proton (938 MeV). 

B. The Weibull Function Fitting Method (WEIB) 

Xapsos, et al. 2 applied a Weibull function to several SPEs that has the mathematical integral form: 


<F = ® 0 * exp (-k*E“) 


where O can be either the proton fluence or proton flux having energy that exceeds a threshold energy E (MeV). The 
constants 4> 0 , k, and a are determined by a nonlinear regression fit to the energy spectrum of interest. The 
differential energy spectrum is easily found by differentiating with respect to E. 

5O/0E = <P 0 *k* a*E“"'*exp (-k*E“) 

They found that the 3 fitting parameters can easily be extracted from either of the above two equations using a 
standard Marquardt-Levenberg nonlinear regression routine 7 ' 8 . For several SPEs analyzed by these authors, they 
show the following Weibull parameters: 
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Table L Weibull Fitting Parameters for Several SPEs 2 


Onset Date 

<P P 

k 

a 

r 

4 August 1972 

2.455x10'° 

0.0236 

1.108 

0.9% 

12 August 1989 

1.622x10" 

1.166 

0.4015 

0.994 

29 Sept 1989 

3.631x10'° 

0.877 

0.3841 

0.998 

19 October 1989 

1.230xl0 12 

2.115 

0.2815 

0.998 

23 March 1991 

1,660x10" 

0.972 

0.4410 

0.993 


In two other papers, Kim, et al. 9 ' 10 analyzed a number of SPEs some of which were GLEs using a Weibull 
distribution of the form 


d> (>E) = P 0 * exp (-a*E b ) 

We compared the results of two of the events in Table 1 above (29 Sept. 1989 and 19 Oct. 1989) with those given in 
ref. 9, Table 4, which is shown below: 


Table 2. Weibull Fitting Parameters for Several SPEs 9 


SPE 

0> (>30 MeV) 

Po 

a 

b 

9/29/1989 

1.35 * 10 9 

1.90 * 10 10 

0.5493 

0.4643 

10/19/1989 

4.23 * 10 9 

1.91 * 10 12 

2.5677 

0.256 


It is also noted that Kim, et al. 9 fitted the SPEs out to a maximum energy value of 1000 MeV (1 GeV). 

C. The Band Function Fitting Method (BAND) 

The fitting method we suggest is that of Tylka and Dietrich 3 and Atwell, et al. 4 where ~70 GLEs dating back to 
1956 were fitted using a Band function 5 , (thoroughly discussed in the above two references) which is a double 
power law in proton rigidity. This method fits actual solar proton data over the entire SPE energy spectrum ranging 
from 10 MeV to 5 GeV. 


III. 12 Ground Level Enhancements (GLEs) 

The 12 GLEs utilized for analysis in this paper are: 

> 1956 February 23 

> 1960 November 12 

> 1960 November 15 

> 1961 July 18 

> 1968 November 18 

> 1971 January 24 

> 1989 September 29 

> 1989 October 19 

> 1997 November 6 

> 2000 July 14 

> 2001 November 4 

> 2003 November 2 

These 12 events were selected, since they were analyzed using a Weibull function by Kim, et al. 9 and provided a 
means of comparison with the Band function. The following plots show a comparison of the three fitting methods. 
Band, exponential and Weibull, for both the integral and differential fitted spectra for the 12 GLEs. Note that we 
also show a comparison of the combined 1960 November 12 and 1960 November 15 GLEs. 

In every case we find that the exponential in proton rigidity fitting method grossly under-estimates both the Band 
and Weibull fitting methods. 
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Comparing the Band method with the Weibull method, we find considerable variation using the Weibull fitting 
method, which we find rather puzzling. The only reasonable agreement was for the 1956 February 23 GLE, and this 
may ju st be a fortuitous fluke. Figures 1. (a) — (m ) compare the integral and differential spectra for thc_GLFs. 
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Figure 1(b) 1960 Nov 12 GLE 
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Figure 1 (c) 1960 Nov 15 GLE 
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Figure 1 (d) 1960 Nov 12-15 
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Figure 1 (e) 1961 July 18 
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Figure 1 (f) 1968 Nov 18 
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Figure 1 (g) 1971 Jan 24 
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Figure 1 (h) 1989 Sept 29 
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Figure 1 (i) 1989 Oct 19 
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Figure 1 (j) 1997 Nov 6 
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Figure 1 (k) 2000 July 14 
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Figure 1 (1) 2001 Nov 4 
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Figure 1 (m) 2003 Nov 2 


IV. Radiation Exposure Calculations 


We now present radiation exposure calculations for the 12 GLEs plus the exposures for the combined 12 & 15 
November 1960 events. Using the NASA LaRC HZETRN 20 10 6 code and the differential spectra for the dozen plus 

1956 Feb 23 GLE - Skin & BFO Exposures 
Band, Exponential & Weibull Spectra 
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Figure 2 (a)”” 1956 Feb 23 GLE 
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GLEs, the skin and BFO exposures were computed and compared with the 30-day and annual NASA crew exposure 
limits. These are shown in Figures 2. (a) - (m). The green dash line indicates the amount of aluminum shielding 
required to not exceed the BFO 30-day limit for some of the events. 



Figure 2 (b) 1960 Nov 12 GLE 


1960 Nov 15 GLE - Skin & BFO Exposures 



Aluminum Shielding, g/cm 2 
Figure 2 (c)^960 Novl5 GLE 
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1960 Nov 12-15 GLEs - Skin & BFO Exposures 
Band, Exponential & Weibull Spectra 



Aluminum Shielding, g/cm 2 
Figure 2 (d) _ 1960 Nov" 12-15 


1961 July 18 GLE - Skin & BFO Exposures 
Band, Exponential & Weibull Spectra 



Aluminum Shielding, g/cm 2 
Figure 2 (e) _ 1961 Julyl8 
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1989 Sept 29 GLE & ESP - Skin & BFO Exposures 



Aluminum Shielding, g/cm 2 
Figure 2(h) 1989 Sept 29 


1989 Oct 19 GLE & ESP - Skin & BFO Exposures 



Figure 2 (i) 1989 Oct 19 
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1997 Nov 6 - Skin & BFO Exposures 
Band, Exponential & Weibull Spectra 
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Figure 2 (j) 1997 Nov 6 
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2000 July 14 GLE & ESP - Skin & BFO Exposures 
Band, Exponential & Weibull Spectra 
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Figure 2 (k) 2000 July 14 
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Figure 2 (1) 2001 Nov 4 


2003 Nov 2 GLE - Skin & BFO Exposures 
Band, Exponential & Weibull Spectra 



Aluminum Shielding, g/cm 2 
Figure 2 (m) 2003 Nov 2 
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V. Conclusions 

A number of conclusions can be made from the previous assessments for the 12 GLEs selected for this study: 

1) The 12 GLEs selected for this study were made since it provided a means of comparison with the Band 
and Weibull fitting methods. 

2) Three spectral fitting methods, Band, exponential and Weibull, were compared; in all 12 cases the 
exponential in proton rigidity fitting method grossly under-estimated the other two methods and 
consequently grossly under-estimated the dose equivalent. 

3) The Band method completely describes the proton spectrum over the entire energy range from 10 MeV 
to 5 GeV. 

4) Only one case, 1956 February 23, showed a reasonable comparison between the Band and Weibull 
spectral fitting methods; for the other 11 cases the Weibull method varied considerably when compared 
with the Band method. 

5) Skin and BFO dose equivalent values for the 12 GLEs were computed using the NASA Langley 
Research Center high energy particle transport/dose code; the values were compared with the NASA 
crew exposure limits for 30-days and an annual limit. 

6) Of the 12 GLEs analyzed, four GLEs (1960 November 12, 1968 November 18, 1971 January 24, and 
1997 November 6) produced an order of magnitude or higher proton fluences, when using the Weibull 
method: 

a. For the 1960 November 12 event, the Weibull method produced ~a factor of 10 higher proton 
fluence for energies <100 MeV when compared with the Band method. 

b. For the 1968 November 18 event, the Weibull method produced ~a factor of 10 higher proton 
fluence for energies <400 MeV when compared with the Band method. 

c. For the 1971 January 24 event, the Weibull method produced ~a factor of 10 higher proton 
fluence for energies <200 MeV when compared with the Band method. 

d. For the 1997 November 6 event, the Weibull method produced ~a factor of 10 higher proton 
fluence for energies <1000 MeV (1 GeV) when compared with the Band method. 

7) Seven of the 12 GLEs required additional aluminum shielding ranging from 2.5-22.0 g/cm 2 aluminum 
to stay at or below the 30-day BFO crew exposure limit as shown in Figs. 2 (a), 2 (b), 2 (d), 2 (h), 2 (i), 
2 (k), and 2 (1). 

8) We strongly recommend that the Band function spectral fitting method be utilized by the 
scientific community for analyzing any future GLEs. 
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